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An Explanation of the Principal Cause of the large Errors at present 
existing between the Positions of the Moon deduced from Hansen’s 
Tables and from Observation; and the Cause of an Apparent 
Increase in the Secular Acceleration in the Moon’s Mean Motion 
required by the Tables , or of an Apparent Change in the Time 
of the Earth’s Potation. By E. J. Stone, M.A.. F.R.S. 

The errors in the Lunar Theory have been traced to the 
effects of changes in the unit of time , or in the point of refer¬ 
ence from which our theoretical longitudes are measured, which 
have, apparently unconsciously, been introduced from time to time 
into astronomy with changes in the adopted data. 

The argument is clearly seen by a consideration of the differ¬ 
ent expressions for the longitudes of what may be called the 
mean Sun which have been adopted for the determination of the 
sidereal times at mean noon. 

If B, H, and Y denote the longitudes of the mean Sun, 
according to Bessel, Hansen, and Le Terrier, we have for 1850 
Jan. 1, Paris mean noon +2 t. 

o / // 

Bs=28o 46 36-12 + 1296027*618184 . £ + 0*0001221805 . t 2 
H = 280 46 43*20+ 1296027*674055 . £ + 0-0001106850 . t 2 
V = 28o 46 43-51 + 1296027*678400 . £ + 0*0001107300 . £ 2 

o a 
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In all these expressions the unit of time has been supposed to be 
a Julian year of 365’25 mean solar days. 

The constant differences, 7 ;/, o8 and 7 //# 39, in B —H and B —Y 
are not unimportant, for they introduce abrupt changes in the 
record of time ; but the differences in the coefficients of t and t 2 
show that the same unit of time cannot have been adopted in 
these expressions. The measure of time must be continuous, and 
the point from which our R. As. are measured, definite. Let, 
therefore, 1 and (1 +a?) be the units of time in B and H. Then 

1296027-618184. 6 + 0*0001221805^ 

== 1296027*674055 . t(l +£’)+ 0*0001106850,. t 2 (l + tf) 2 . 


If, therefore, 

n' = 1296027*674055 

o-o5587i , Q-OOOII4955 

n' w 

To reconcile B and H, therefore, x must contain a variable term. 
Similar remarks apply to the difference between B and Y. 

Now let N be the Moon’s mean motion referred to 1 as the 
unit of time, and (N -f £N) the Moon’s mean motion referred to 
(i+as) as the unit of time. Then 

(R + 5N)(i+ar) = N, 

and 

*.5N = X{0-055871 00000114955 ■< J }=o"747 • • 

But Hansen determined his mean motion of the Moon so as 
to force an agreement between his theory and observations 
reduced with Bessel’s unit 1, and his Tables therefore repre¬ 
sented the observations well for many years whilst 1 was adopted 
as the unit of time; but directly the unit of time was changed 
by the adoption either of H or Y, then the effects of the 
erroneous determination of the Moon’s mean motion by Hansen 
become apparent. 

The effect of the change of unit is also shown in the com¬ 
parison of Le Yerrier’s Solar Tables with observation, but of 
course only to about the thirteenth part of the amount shown by 
the Lunar Tables. 

Le Yerrier’s Tables have been employed in the computation 
of the Sun’s positions in the Nautical Almanac since 1864; the 
Greenwich observations available to me only extend to 1880. 
The following table exhibits the mean excess in longitude of Le 
Yerrier’s tabular places over those deduced from observation. 
It will be seen that the errors are clearly progressive, and the 
change per annum appears roughly to be about +o /7 *o52, a 
result confirming, I think, my opinion of the undesirability of 
the increase of the mean motion of the Sun by +o //, o63. 
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Year. 

in 

Longitude. 

No. of OK'. 

Year. 

in 

Longitude. 

No. of ' >1 is. 

1864 

4- 0*025 

122 

1873 

it 

+ 0*408 

124 

1865 

+ 0*136 

222 

1S-4 

+ 0*762 

106 

1866 

— 0*109 

109 

1S75 

+ 0*728 

no 

1867 

— 0-203 

93 

1876 

+ 0*675 

112 

1868 

+ 0245 

131 

1877 

+ 0*683 

87 

1869 

+ 0*272 

98 

1S7S 

+ 0*816 

85 

1870 

+ 0*122 

126 

1S79 

+ 0*748 

77 

1871 

1872 

+ 0*306 

-OO 37 

113 

120 

1SS0 

+ 0861 

“5 


‘Table showing the agreement of Hons- n + Tr viTi actual observation before the 
change in the unit of time . and the /a ~ gr:<-'iee >. rear after the change in the 
unit of time introduced in 1864. 


Year. 

Mean Value 
for Year 
(H — 0). 

No. of Obs. 

Y. or. 

Mean Value 
for Year 
(H — 0 ). 

No. of Ob; 

1847 

n 

+ 1*066 

107 

*^55 

// 

-O-871 

94 

1848 

+ 0-195 

113 

1S56 

-0-956 

93 

1849 

-0*465 

109 

1857 

-1-858 

"5 

1850 

— 0*282 

114 

1S5S 

-i -975 

130 

1851 

-1*295 

109 

1S60 

not compared 


1852 

-0*917 

109 

1S61 

? 5 


1853 

-1*629 

105 

1S62 

— 2-829 

§5 

1854 

—1*684 

121 

1S63 

1*606 

ir 5 

Here change in 

the unit of time 

took place. 


1864 

+ 0*121 

hi 

1S73 

+ 8*239 

94 

1865 

+ 1*271 

108 

1S74 

+ 9-294 

1 17 

1866 

+ 2*142 

107 

1 —75 

+ 9*867 

94 

1867 

+ 3*480 

89 

1S70 

- 9-755 

82 

1868 

+ 4*117 

125 

1S77 

+ 9-234 

105 

1869 

+ 4'277 

94 

1S7S 

4-8*219 1 

90 

1870 

+ 4-828 

102 

1S79 

+ 9*631 

§5 

1871 

+ 6-955 

106 

18S0 

+ 10*265 

IOI 

1872 

+ 7*309 

109 





G G 2 


© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 


Downloaded from http://mnras.oxfordjournals.org/ by guest on June 9, 2015 





1883MNRAS..43..335S 


338 Mr. Stone , An Explanation xliii, 7, 

To get rid of the effects of the diurnal rotation I have re¬ 
ferred ever jthing to the meridian of some Observatory, as Paris, 
and discussed the changes of sidereal times at mean Paris noons ; 
but as there is apparently some difficulty in following this reason¬ 
ing, which is based upon the principle of continuity, I give a 
direct proof of the legitimacy of the method adopted in throwing 
the changes in the law fixing the position of the point from 
which our theoretical R.As. are measured, relatively to the Paris 
meridian at mean noon, on the time in the case of the mean Sun, 
from whose motion our ideas of mean solar time are derived,, 
and on the determinations of the mean motions for the Moon 
and the planets. 



If we assume the time of the Earth’s rotation on its axis to 
be constant—an assumption to be tested by the results—we can 
always determine the position, relatively to our meridian, of a 
fixed point on the equator at any definite time. Let this point 
be 0 . This is a simplification merely to get rid of the complica¬ 
tion of ideas arising from the diurnal motion. 

Let m be the mean Sun by whose motion our time is to be 
measured, and at time t, expressed in any definite unit, let its 
position be at m, so that 

0 m « nt 


(ji is a definite numerical quantity, since the rotation of the Earth 
on its axis is supposed to be assigned) ; and at the same time 
let P be a planet, and 

OP=A + n't. 

For practical convenience we measure our angles from y, where 

Oy = C -{-at + st 2 , 
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0, a, and s being assigned constants. 

7 m - C -f (a -f «}* + s£ 2 
7P« C 4- A -f (s + + st 2 . 

But suppose, instead of referring our angles to y, we refer them, 
in our theoretical investigations, to y n such that 

0 % = C + 5 C - ( ? - 0 a): 4- (* 4- 5 5) ; 

then 

( 1 ) 7 ,m = C + 8 C *f (a -r ;i -f 5 a)/ + (s + 8 s)£ 2 , 

t having the same value as before. 

(2) 7 y P = C + 5C + A-r(f/*f3c-f ;i')A + (s + 8 s)£ 2 . 

Now, if we know of the changes SC, Ss, and allow for 
them, we shall correctly determine the relative motion of P 
measured from y ,; but if we do not know of these changes, and 
do not allow for them, how shall we force an agreement between 
observation and theory ? Clearly we shall have to throw the 
error on the time measured from (1), and on the mean motions 
of the planets in (2), and we shall thus partially force an agree¬ 
ment, but not a complete one, between observations referred to y 
and theory referred to y r 

First, let t be the true time deduced either from ym or y; 
and let (^ + $0 be the observed time deduced from (1) when 3 C, 
da, and ds are neglected, or when y J m is forced to agree with 
observations referred to y; and let B be the epoch, x the mean 
motion, and yt 2 the empirical secular correction required to make 
y / P, deduced with the time (t+dt), agree with yP; we then 
have 

7 /P~b4 <x{t + 8<Q + (s + 8s-f y)t~ 

— 0 + A. + (jet + + st~. 


But from (1) we have 


and 


0 = 3 C + (a + n)8 t + oc(j . r -f os. i~ 
g ^ — (5C 4 - 8 a . £ 4 - 8s. i ~j 

a + n 


-r, 0 ^.5C , f 8 a \ { 

7,P = B - bxl I-- \t 4 ( s 

' a + n \ a -\-n) \ 

~ C + A 4 - [a 4- ?i f ) . t 4 - st 2 

B—A 19 = C + A 

a 4 - n 

x(i — JPP \ = (a + n f ) 

\ a + n/ 


y + 8s 


a +n) 
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, * 5s. r 

fi+^ + 05 -- s 

« + n 



8 a \ 
a-r7i) 


B=C+A+5C 

\a + n J 


y + 8 s = 8 s 



The correction required to y P deduced directly from the 
erroneous time to make it agree with the observations referred to 
y is therefore 

-(5C + 5a.* + $s A 2 ) ( c Et^\, 

\a + # ) 

which is the condition upon which the investigation in my paper 
is based as deduced directly from the consideration of the 
differences in the units of time necessary to secure continuity 
in the adopted expressions for sidereal time at mean noon. The 
possibility of such errors arises from the indirect way in which 
our theories are compared with observation. Numerical values 
of the constants which fix the relative changes in the planes and 
points of reference and the mean motion of Sun have to be 
adopted and carried into the theories ; but the observed R.As. 
are deduced from the times of transit by clocks whose errors are 
determined by the differences between the times of transit by the 
clock and R.As. measured from the point of intersection of 
the ecliptic and equator, which point must, for accurate work, 
be the same as that adopted in the theoretical expressions. I 
have shown, however, by the discordant values of B, H and V, 
that no definite point, much less the true point, has been practi¬ 
cally adopted in our theories. 


z 
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The breaks of continuity in our measures which have neces¬ 
sarily taken place with the adoption of different units through 
the introduction of small chancres in the Sun’s mean notion, and 
in the secular terms which ex press the relative motions of the 
planes of reference, are clearly shown as follows :— 

Let x, y, Z be the projections of three rectangular axes fixed 
in space, xy being the plane of the ecliptic and x the point of 
intersection of the ecliptic and equator at the instant from which 
onr time is measured ; and at the time t let y 0 AB be the plane of 
the equator, yN the plane of the ecliptic, C the projection of the 
polar axis of the Earth, and CA the meridian of the Paris Obser- 
tory. 

Let 

ZC = 0; 7 c A « C ; 7 0 A = <j>; 

n = angular velocity of th-_- Karrh on its axis. 


Then 


And 



This equation is exact. 


f = G + nt + 



e 


where C is some constant to be determined from observation. 
Prom the theory of the motion of the Earth about its centre of 
gravity, we have, to the second power of /. 

(1) $ = at + bt‘- V. 

(2) e = a 0 +/ C - C .: 

where a, b, w 0 , / are constants, th and O periodic functions. 

. cos 9 .dt —at cos (o Q + cos v c + ¥ . cos a Q , 
clt 

and 

<p = C + {71 + a cos « 0 ) . t - C- cos oo c -W . cos a> 0 . 

But from the relative motion of the Sun around the Earth, we 
have 

(3) — fit; (4) t+et 2 , 

where a, ( 3 , y, and c are constants. 

Prom (i), (2), (3), and (4) we obtain 

y Q y = ut 4 - u't 2 ; where u and vJ are constants. 
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But 

7A = <£ — 7 c 7 = s, suppose; 

5 = C + (?£ + a cos w 0 — w)jf + (b cos w 0 — u') t 2 + W cos » 0 . 

0 

This equation fixes the law of variation of y the point from 
which the R.As. are measured, relatively to A, a fixed point on 
the Earth’s surface ; and nothing but confusion can result from 
any alteration in this expression unless the changes in the con¬ 
stants by which it is produced are consistently carried through 
all our work. 

Aow suppose the mean longitude of the Sun 

— L 0 + .NY + A t~. 

Let 

w + (2 COS W 0 365*25 X27r; 


this fixes the unit of time where 


11 = angular Telocity of Earth on its axis, 
a = luni-solar precession, 

« 0 = obliquity at 1850, Jan. I, 
u — intercept on Equator between ecliptic-planes, 
n' = Sun's sidereal mean motion in longitude, 
p = general precession. 

N = n' +p. 


All these quantities, when expressed numerically, have to be 
referred to a Julian year of 365*25 mean solar days. But since 
the rotation of the Earth on its axis, and the Sun’s mean motion 
in longitude are physical facts, we must have 



where C is a constant to be determined from observation; and 
this relationship must hold good whatever unit of time he adopted , 
and be carried into the equation of condition which fixes the unit 
of time from the relative motions of the Earth on its axis, and 
the Sun’s mean motion in It.A. 

Similarly we may write 

a cos & 0 — n — p — Ant; 

where A is some constant. 

Our equation of condition which defines the unit of time 
then becomes 

»'(C—I + A) = 365*25 x 2 ir. (2) 

This equation shows at once that we cannot change n r without 
changing onr unit of time. We can of course write 

.(»' + (0 -1 + A) = 365-25 x 2 tt, (3) 
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and still call the new nnit of time a Julian year , bnt in this case 
the length of a Julian year in (3 ) is not the same as in (2); and 
if t l and t 2 are the lengths of the two so-called Julian years, we 
must have 

50 ; ( 4 ) 


therefore 





(5) 


From which it directly follows that if a period be expressed by 
t Julian years under system (2), it will exceed the period 
expressed by t Julian years under system (3) by 

I now proceed to show how it has been frequently assumed 
that the error due to a change bn* is only 

_ 8n' ^ . ... .! f 5/i r . 

n " ** n- 


Let n and n + on be the angular velocities of the Earth on 
its axis, corresponding to the two values h and (V + In') of 
the Sun’s mean motion. Then 


n = 365*25 x 2 7 r + n f +p + ii—a e s a-. ; (7) 

but it is then usually assumed that 

n + Bn =365*25 x 27 t + n r + + u ~ a cos (8) 


from which of course would follow 


Bn — 5 '/ 


and 

Bn _ B rt _ 
n n 

bnt upon this I must remark that since r, and t 2 cannot both be 
put —1, we must have instead of (8), the equation of condition 

m + Bn — 365-25 x 27 t(i 4- #) + »' + Bn' 4- (1 -a) (p + u — a cos co Q ). (9) 


v(The change of unit in n f is allowed for in bn'.) 

.*. Bn — 365*25 x 2ttx + 5 n e -b x (p -f 11— a cos w Q ) ; 
5 n — “-(365*25 x 27 r + n’ —p -f u — a cos cc Q ) ; 



. Bn _ Bn' 

• * ~~ 7} 

n n 
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or tlie change of unit of time is as 


i : i + 


8 n f 


l 1 as before proved. 

Since the epoch from which t is measured is not yet fixed, let 
it be so chosen that 


Then 


C = L n . 


s = L 0 + N£ + 365-25 x 27 r x t+ (b cos 03 o —u')t 2 + W .cos o> 0 . 


If the epoch'thus chosen be called 1850, Jan. 1, Paris mean 
noon, we have for the computation of the distance y from A at 
the different Paris mean noons, 

S Q = L 0 + Nt + (b . cos a> 0 — 'ii') .t 2 + ¥. cos co 0 . 


The non-periodic part of S 0 only differs from the mean longitude 
of the Sun by 

{X — (b. cos eo 0 — u')} . t 2 . 

If, therefore, we call 

L=L 0 + -f (b . cos co 0 — u f ) . t 2 


the longitude of the mean Sun, the mean longitude of the 
Sun 


= L + {A — (b . cos . t 2 . 


Here 


| A — (5 cos cc a — u ')}. t 2 


is a small secular term which, with the data adopted in Le 
Yerrier’s Tables— 


= — 0-00003113 . t 2 = —0*311 



If it should be found absolutely necessary, in the progress of 
science, to introduce corrections SJST and S\ to N and X in order 
to represent the solar observations, we must keep 


L = L 0 + ISiY-f (6 cos <a 0 —u’)t 2 


still unchanged, but put the mean longitude of the Sun 


= L + 8M + ( 5 a-o*oooo 3 H 3 ) .t 2 . 


If the constants determining the motions of the planes of 
reference and the mean motion of the Sun are accurately known, 
we shall, in this way, have for a very long period the point of 
reference y coincident with the point of intersection of the 
ecliptic and equator; but even if this condition is not quite 
satisfied, we must not change y by changing the law of variation 
of Ay until we deliberately change our unit of time, and carry 
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logically through all our work the consequences of that change. 

If the difference 0*3113 between the longitude of the 

mean Sun and the mean longitude of the Sun be neglected (and 
this has usually been done),it would lead to the apparent necessity 
of an increase of about 4"-16 in the secular acceleration of the 
Moon’s mean motion. This quantity agrees very closely with 
the supposed discordance between theory and observation. It 
would appear, therefore, that the discordance between theory 
and observation which has recently been thrown on the Earth’s 
rotation no longer exists to any proved sensible extent, and, so 
far as we know at present, the time of the Earth’s rotation is 
constant. 

[This abstract was drawn up for the Society’s publications 
by request, after a verbal explanation of the paper, which was 
given at the meeting of April 13. 1SS3." 


On the Computation of the Eccentric Anomaly , Equation of the 
Centre and Radius Vector of a Pin net, in Terms of the Mean 
Anomaly and Eccentricity. By J. Morrison, M.D., M.A., 
Assistant on the American Ephemeris and Nautical Almanac , 
Washington, D.C. 

(Communicated by the Foreign Secretary.) 


Several methods of computing the eccentric anomaly of a 
planet from the eccentricity and mean anomaly have already 
been given, in the Monthly Notices and elsewhere, by many very 
able mathematicians. Some of these methods involve consi¬ 
derable labour even after the formulas have been deduced, while 
others give only rough approximations. The following method 
of treating the subject is due, I believe, to the late Professor 
Hansen, as indicated by him in the Ahhandlungen der Sdchs- 
ischen Cesellschaft der Wissenschaften, Band II. It gives the 
eccentric anomaly with great accuracy and facility by means of 
rapidly converging series, and in most cases with very little 
labour, especially after the coefficients of the several terms of the 
series have been computed for each planet. I purpose in this 
short paper to develop the subject more fully than is done in the 
work referred to, and to apply the results to the computation of 
the eccentric anomaly in the case of each of the primary planets, 
as well as of the equation of the centre and the radius vector. 

Let 


and 


M = the Mean Anomaly, 

E = the Eccentric Anomaly, 
v — the True Anomaly, 
e — the Eccentricity, 

e = the Napierian base- 
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